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Abstract 

Self-scaled barrier functions on self-scaled cones were introduced through a set of axioms in 
1994 by Y. E. Nostcrov and M. J. Todd as a tool for the construction of long step interior point 
algorithms. This paper provides firm foundation for these objects by exhibiting their symmetry 
properties, their intimate ties with the symmetry groups of their domains of definition, and 
subsequently their decomposition into irreducible parts and algebraic classification theory. 
In a first part we recall the characterisation of the family of self-scaled cones as the set of 
symmetric cones and develop a primal-dual symmetric viewpoint on self-scaled barriers, results 
that were first discovered by the second author. We then show in a short, simple proof that any 
pointed, convex cone decomposes into a direct sum of irreducible components in a unique way, 
a result which can also be of independent interest. We then show that any self-scaled barrier 
function decomposes in an essentially unique way into a direct sum of self-scaled barriers 
defined on the irreducible components of the underlying symmetric cone. Finally, we present 
a complete algebraic classification of self-scaled barrier functions using the correspondence 
between symmetric cones and Euclidean Jordan algebras. 
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1 Introduction 



In recent years a theory of interior-point methods for hnear, semidefinite, and second-order cone 
programming has been developed within the unified framework of self-scaled conic programming. 
The origins of this theory can be traced to the works Q . The importance of the problems 

which can be cast in this framework, and the fact that it is possible to develop efficient long-step 
interior-point methods for these problems have contributed to the popularity of the subject in the 
optimization community and beyond. 

In order to facilitate our exposition, we consider the following pair of convex programs in conic 
duality 

(P) infix, So) (D) mi(xo,s) (1.1) 

xe{L + xo)nK s e {L^ + So) n K*. 

Here i5 is a finite dimensional Euclidean space equipped with an inner product (•, •), L is a linear 
subspace of E, and its orthogonal complement. The cone K is a. regular (closed, convex, 
pointed, solid) cone, Xo G int(i^), and Sq S int{K*), where K* is the dual cone 

K* -.^{seE: {x, s) > 0, Va; e K}. (1.2) 



Interior-point algorithms can be used to solve (1.1) over any regular cone, provided one has 
a self-concordant barrier function F(x) defined over the interior int(_R') of K. The reader may 
consult the authoritative monographs |l^, |l^ for a detailed treatment of self-concordant functions 
and interior-point methods. In a generic self-concordant barrier function, one has control over the 
behaviour of the Hessians F"{y) only when y lies in the local ball {y : {F"{x){y — x),y — x) < l}, 
leading to short-step interior-point methods. Although these methods have a polynomial running- 
time guarantee, they tend to be less efficient linear programming solvers in practice than long-step 
interior-point methods. The theoretical basis for this latter type of algorithm is the fact that the 
self-concordant barrier function 

n 

X — In Xi 
1=1 

has additional properties which make it possible to control it in all of int(_ftr). 

In iQ, Nesterov and Todd isolate two properties of the barrier — X)"=i ^^^i which they identify 
as being responsible for makin g th e lon g-st ep approach succeed in linear programming. They 



generalise these properties (see {lA) and (L5)) and call self-concordant barrier functions satisfying 
these conditions self-scaled. Since these properties also impose certain conditions on the domain of 
definition of such functions, Nesterov and Todd call the closures of such domains self-scaled cones. 
For convenience, we recall these concepts here: 

Definition 1.1. Let _ftr C £' be a regular cone. A self-concordant barrier function F : int{K) M 
is called self-scaled if F"{x) is non-singular for every x G K, F is logarithmically homogeneous, 
that is, there exists a constant v > such that 

F{tx) = F{x) - i^int, Vx G int(i^), t > 0, (1.3) 

and if F satisfies the following two properties 

F"{w)x G int{K*), Vx, w G int{K), (1.4) 
F^ {F"{w)x) = F{x) - 2F{w) - v, Vx, w G int(i4:). (1.5) 

If K allows such a barrier function, then K is called a self-scaled cone. 



2 



The dual barrier : int(_fi') — > R that appears in the last Axiom (1.5) is defined as F,,{s 



sup{ — (a;, s) — F{x) : x e int(_ft')}. Theorem 3.1 in jlj] states that (1^) can be strengthened to 



Theorem 1.2. If x ^ int(A') and y G int(iir*), then there exists a unique point w G int(-R') such 
that 

F"{w)x = y. 

Moreover, if w € int{K) then F" {w){K) ^ K*. □ 

The point w is called the scaling point of x and y. The last statement is a consequence of the 
first part of the theorem and of Equation (3.2) from Nesterov and Todd's paper ||lj]. We reproduce 
this formula here for convenience: 

F"{x) = F"iw)Fi'{F"{w)x)F"{w). (1.6) 



See also Lemma 2.3 below, where this formula reappears and where we give a proof of this important 
identity. 

We would like to mention that Rothaus using rather elementary tools, proves a number 
of results which are useful in Section of this paper. Two key results are Q, Theorem 3.12 



and Corollary 3.15, p. 205. These results imply Theorem 1.2 for the universal barrier function, a 



special self-concordant barrier function defined by Nesterov and Nemirovskii [tl3l which is further 



discussed below. Theorem 1.2 is an independently discovered extension of Rothaus's result to all 



self-scaled functions. It can be shown that, when suitably modified, all results of Section III in 



|18| can be extended to general self-scaled barriers. In particular this is true for Theorem 1.2 



Nesterov and Todd ||l^, 15 demonstrate that self-scaled barrier functions can indeed be used 
to develop various long-step interior-point methods for linear optimization over self-scaled cones, 
in particular for semidefinite programming and for convex quadratic programming with convex 
quadratic constraints. 

Inspired by the paper of Vinberg O. Giiler H develops the relationship between the uni- 
versal barrier function of Nesterov and Nemirovskii |13] and the characteristic function of the cone 

^k{x) := / e-<^^^>dy, (1.7) 

J K* 

introduced in 1957 by Koecher, see |l^. In particular, Giiler |^ shows that if if is a regular cone, 
then the universal barrier function U [x) for K satisfies the equation 

U{x)=ci\nipK{x)+C2, (1.8) 

for some constants ci > and C2. 

Through Giiler's paper Q the concepts of homogeneous cones, homogeneous self-dual cones 
(or symmetric cones), Euclidean Jordan algebras and Siegel domains as well as the classification 
theory of symmetric cones and Euclidean Jordan algebras, known to mathematicians since 1960 
and 1934 respectively, were first introduced into the in the interior-point literature. The interested 
reader can find a complete treatment of these classification results in the book of J. Faraut and 
A. Koranyi See also (lo) for a different treatment of some of the same topics. Because of their 
importance for this paper, we recall some of the concepts mentioned above. 

Definition 1.3. Let K Q E he a. regular cone. The automorphism group of K is the set of all 
non-singular linear maps A : E ~f E that leave K invariant, i.e., 

Aut{K) := {A e GL{E) : A{K) = K} 
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The cone K is called homogeneous if Aut(_R') acts transitively on if, that is, given arbitrary points 
x,y E K, there exists a map A E Aut(K) such that Ax — y. The cone K is called self-dual if E 
can be endowed wit h an inner product such that K* = K where K* is defined with respect to this 
inner product, see (L2). The cone K is called homogeneous self-dual if K is both homogeneous 
and self-dual. 



Homogeneous self-dual cones are also called symmetric cones in a terminology which we 
shall adopt in this paper. 

The motivation behind |l^ in contrast to is rather different: While the first two papers 
deal with long-step interior-point methods and regular cones on which such methods can be 
designed, the latter one deals with the universal barrier function and the symmetry properties 
of regular cones, both in the group theoretic and the duality theoretic sense. Shortly after the 
announcement of the paper Giiler discovers that the families of self-scaled cones and of 
symmetric cones are identical, thus establishing a connection between these two previously distinct 
ideas. 



As mentioned earlier, symmetric cones are fully classified in the theory of Euclidean Jordan 
algebras, see Q and the references therein. According to this theory, each symmetric cone 
has a unique decomposition into a direct sum of elementary building blocks, so-called irreducible 
symmetric cones, of which there exist only five types. Three examples of symmetric cones are of 
particular interest to the optimization community: The non-negative orthant K = R" , the cone 
K — of n X n symmetric, positive semidefinite matrices over the real numbers, and the Lorentz 
cone K = {(|^) E R"+^ : r > ||a;||2}. The general self-scaled conic optimization problems associ- 
ated with these cones are respectively linear programming, semidefinite programming and second 
order cone programming. The latter can be seen as a reformulation of convex quadratic program- 
ming with convex quadratic constraints. Considering more general symmetric cones, one can treat 
linear optimization problems with mixed linear, semidefinite and convex quadratic constraints in 
a single unified framework. 

Motivated by Q and by the fact that only a small number of examples of self-scaled barrier 
functions are explicitly known, Hauser develops a partial algebraic classification theory for self- 
scaled barrier functions in a chapter of his thesis |Q and later announces a report based on 
these result. Hauser shows that any self-scaled barrier over a symmetric cone K has an essentially 
unique decomposition into a direct sum of self-scaled barriers defined on the irreducible summands 
of K. Using this decomposition, he classifies in particular the family of isotropic self-scaled barrier 
functions which are characterised by rotational invariance. The insight gained from a lemma leads 
Hauser to conjecture that all self-scaled barrier functions defined on irreducible symmetric cones 
must be isotropic. Hauser points out that this conjecture, if true, would settle the classification 
problem of general self-scaled barrier functions. This conjecture also implies that all self-scaled 
barriers over irreducible symmetric cones are of the form ci In ipx + C2, where ci > and C2 are 
constants. It follows from this theory that the full set of self-scaled barrier functions is readily 
known, and that all of these functions are just minor transformations of the universal barrier 
function. 

In a second report |9|, Hauser proves this conjecture in the special case where K is the positive 
semidefinite cone, i.e., he shows that all self-scaled barrier functions for use in semidefinite program- 
ming are isotropic and essentially identical to the standard logarithmic barrier function. Starting 
from first principles, Hauser shows that the orientation preserving part of the automorphism group 
of the positive semidefinite symmetric cone is generated by the Hessians of an arbitrary self-scaled 
barrier function defined on this cone, see Corollary 4.3. Hauser's solution of the isotropy con- 
jecture in this special case relies on Proposition 3.3, which also forms the key mechanism in the 
proof of Corollary 4.3 in |^. 
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Hauser's Corollary 4.3 is essentially a rediscovery of a result by Koecher, Theorem 4.9 (b), pp. 
88-89 jl^, for the special case of the cone of positive semidefinite symmetric matrices. Shortly 
after the announcement of Hauser's report |^ in March of 2000, Y. Lim ijlll settles the general case 
of the isotropy conjecture by generalising Hauser's proof while refereeing the paper. Subsequently, 
both O. Giiler and S. Schmieta independently of each other and independently of Y. Lim 
prove the isotropy conjecture in the general irreducible case. It is interesting to note that all three 
approaches to the general case of the conjecture, as well as Hauser's approach to the special case 
relevant to semidefinite programming, rely on the same deeper principle provided by Koecher's 
Theorem 4.9 (b) cited above. See also Remark of this article. 

The present article is a revision of Hauser's original report which is based on parts of his thesis 
0, while incorporating the solution to the general problem using Giiler's approach. Schmieta's 



report |1S[] constitutes the first document where a proof of the general classification Theorem 5.5 
became publicly available. 

The rest of the paper is organised as follows. In Section || we reconsider self-scaled cones and 
self-scaled barriers from a symmetric point of view. Section ^ is devoted to certain properties 
of self-scaled barriers which link self-scaled barriers to the symmetry group of their domain of 
definition. These results are needed in later sections. In Section ^, we show that any pointed, 
convex cone has a unique decomposition as a direct sum of irreducible components. This result, of 
which we manage to locate only technically more involved generalisations, may be of independent 
interest. We therefore include a simple proof. We then use this decomposition result to show 
that any self-scaled barrier defined on a symmetric cone K decomposes in an essentially unique 
way into a direct sum of self-scaled barriers defined on the irreducible components of K, which 
also shows that the irreducible components are symmetric cones themselves. This decomposition 
reduces the problem of classifying self-scaled barriers t o th e case where the domain of definition is 
irreducible, a problem we solve in Section |^. Theorem 5^ constitutes the main and final result of 



this paper. We thus present all the essential elements of the theory of self-scaled barrier functions 
in a single document. 

The following basic properties of z^-self-concordant logarithmically homogeneous barrier func- 
tions and their duals will be used frequently in later sections. These properties are easy conse- 
quences of logarithmic homogeneity, see or : 

Proposition 1.4. Let F be a u- self- concordant logarithmically homogeneous harrier function on 
the regular cone K G E, and let x G 'nit{K), s € mt{K*). Then 

i) - F'{x) = F"{x)x e mt{K*), ii) - Fl{-F'{x)) = x, 

lii) F'J{~F'{x)) = F"(a;)-i, iv) {x, -F'{x)) = 

v) F^^\tx) ^t-^F^^\x), Vt>0,fc = l,2, vi) F^{-F'{x)) ^ -V - F{x), 
vii) F{x) + F^{s) > ~v — vlogv — v\og{x, s), 

where F^^\x) = F'{x), F^'^\x) = F"{x) in v). □ 



2 A Symmetric View on Self— Scaledness 

In this section we undertake a study of self-scaled cones and barrier functions while emphasising 
their symmetry properties in a duality-theoretic sense. 

Let F be a self-scaled barrier function on a regular cone X in a finite dimensional Euclidean 
space E equipped with an inner product (•, •). With a given arbitrary point e € 'mt{K) we associate 
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an inner product 

Me {F"{e)u,v). 

The following result is due to Giiler 

Theorem 2.1. The cone K is symmetric, and F is self-scaled under (•, •)e. 
Proof. We have 

K: := {y : {x,y), > 0, Vx G if } = {2/ : (F"(e)x,2/) > 0, Vx G if} 
= {y : (z,y> > 0, Vz G K*} - (if*)* = if, 
where the third equality follows from Theorem |1.2| . Note that 

{F"{x)u,v) = = v)e = {F" {e)Fi' {x)u,v) 

yields F"(a;) = F"{e)F'J{x), or 

= (2.1) 
Theorem O implies that F^'{x){K) = F" (e)-^ F" {x){K) = F"{e)-^{K*) = if, so that F^'ix) G 



Aut(if). Theorem 1.2 also shows that, given any two points u,w G int(if), we can find a (unique) 
point z G if such that F"{z)u = F"{e)v G if*. Therefore, F^'(x){u) = which shows that the set 
of linear operators {F^'{x) : x G int(if )} acts transitively on if. Hence, if is a symmetric cone. 

For the second assertion, note that if s G if* = if, then 

(Fe)*(s) sup {~{x, s)e ~ F{x)} = sup {-(x, F"(e)s) - F{x)} = F,[F"[e)s). 
xeK xeK 

For a;, z G int(if ), we thus have 

iF,)4F:^'{z)x) - F4F"{e)Fi'iz)x) = F,{F" {z)x) - F{x) - 2i^(z) - 



where the second and last equalities follow from (2.1) and (|l.5|), respectively. Consequently, F is 



self-scaled under (•, Oe- □ 



Remark 2.2. Theorem 2.1 shows that from here on we may assume without loss of generality 



that if is a symmetric cone and that there exists a (unique) point e G int(if ) such that F"{e) — I. 



Together with Equation (g_^) this implies that 

F,(.t) = F{x) - 2F{e) -v = F{x) + const, (2.2) 
and invoking ( |l.5| ) once more this implies the identity 

F{F"{w)x)=F{x)^2F{w) + 2F{e), Vx, w G int (if ). (2.3) 
Note that ( ^.3| ) is a criterion that involves only the primal barrier F. Indeed, t his id entity allows one 



to characterise self-scaled barrier functions without invoking Fi, , see Lemma 2^ below. Changing 
a barrier function by an additive constant is of no real consequence, as interior-point methods rely 
on gradient and Hessian information. T heref ore, we can assume without loss of generality that 
F{e) = 0. For the same reason. Equation ( |2.2| ) makes it possible to think of F and as the same 
function. Hence, we no longer need to distinguish between primal and dual quantities - between 
F and i^*, the primal and dual scaling points and so forth. 

We next prove a property of the Hessian F"{w) which will become an essential tool for our 
classification of self-scaled barriers. For all y G int(if ) let us define 

P(y) -.^F'^yy^ 



6 



Lemma 2.3. For all x,w ^ int{K) it is true that 

P{P{w)x) = P{w)P{x)P{w). 



(2.4) 



Proof. Let us define z = P{w)x. Equation ( |2.3| ) implies that for any h £ E we have F{F" {■w){z + 
th)) = F{z + th) — 2F{w) + 2F{e). Expanding both sides of this equation and comparing the t'^ 
terms one gets 

D^F{F"{w)z)[F"{w)h, F"{w)h] = D'^F{z)[h, h], 

or {F"{x)F"{w)h,F"{w)h) ^ (F"{z)h,h). Thus, F" {w)F" {x)F" (w) = F"{z) = F"{P{w)x), and 
(O) follows. □ 



In the proof above, we need only the weaker condition F{F" {w)x) — F{x) + c{w) where c{w) is 
a constant dependent on w. (However, Lemma |2j| below shows that this is equivalent to (|2.3D.) 



Equation (2.4) is a symmetric version of a Formula (3.2) from jlj], see also Equation (|L6) above. 
In accordance with the established tradition in the theory of Jordan algebras we call (2.4) the 
fundamental formula. 

Remark 2.4. We do not have a Jordan algebra connected to F yet, but the fundamental formula 
leads one to suspect that there might be one. In spring 2000, inspired by the work of Petersson 
p6t , Giiler proves that this is indeed the case. Subsequently, S. Schmieta independently 
discovers the same result, following essentially the same steps. Schmieta uses this fact as an 
essential tool to classify self-scaled barriers. As it turns out, the Jordan algebra connected to F 
is already discovered by McCrimmon in his thesis , even without the assumption of convexity 
for F. His proof in turn is a generalisation of Koecher's ideas on uj -domains. Reading both 
works is instructive in delineating the role of convexity. 



The following result provides an alternative definition of self-scaled barrier functions. 

Lemma 2.5. Let K be a regular, self-dual cone. A logarithmically homogeneous self-concordant 
barrier function F on mt{K) is self-scaled if and only if 

F"{w)x e iiA{K), Vx, w e int(i4:), (2.5) 
F{F"{w)x)^F{x)+c{w), Va;,w G int(i^), (2.6) 

where c{w) is a constant that depends on w. 



Proof. Since K is self-dual. Equation (2.5) is equivalent to Axiom (1.4). If F is self-scaled, 
then Equation (2.6) follows from (p.3[). Conversely, assume that F satisfies Equation (2.6). Let 



X, s Cz int(_ftr) be arbitrary points. 



We claim that there exists a point w G int(_ftr) such that F"{w)x = s. Towards proving 
the claim, we consider the optimization problem min{(z*,a;) : (z,s) = 1}, where z* = —F'{z). 
It is well known that the feasible region is bounded, see pj. Corollary 1.1.6, p. 4. We have 
F{x) + F*(z*) > —V — vlogv — v\og{z* ,x) (see Proposition 1.4 vii))^ and F{z) + F*(z*) = —v 
(see Proposition 1.4 wi)), which imply F{x) — F{z) > — t^log:/ — i/log{z*,x). These imply that 
the objective function of the optimization problem goes to infinity as z approaches the boundary 
of the feasible region, and thus the optimization problem has a minimizer z G mt{K) satisfying 
F"{z)x — As for some scalar A. Since F"{z)x,s G int(i^), we have A > 0. The point w = ^/Xz 
satisfies F"{w)x — s (see Proposition 1.4 w)), and proves the claim. 



Next, we claim that 



c{w) = ~2F{w)+2F{e). 



(2.7) 
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Let u £ mt{K) be a point satisfying F" {u)w — e. The fundamental formula (2.4) is a consequence 
of dU) and gives F"{u)P{w)F"{u) = I, or equivalently F"{w) = F"{u)^. From (|U), we obtain 

F(e) + c{w) = F{F"{w)e) = F{F"{ufe) = F(e) + 2c(w), 



or c(w) = 2c(u). Equation (2.6) also implies that 



F(e) = F{F"{u)w) = F{w) + c(u) = F(u;) + -c(w), 
hence proving the claim. 

Using logarithmic homogeneity alone one can prove that F^(w*) = —v — F(w) where w* :— 
—F'(yj) (see Proposition 1.4 wi)). Proposition 1.4 ii) shows that the mapping w i-^ w* is involutive, 
that is, w** = w. These imply F*(w) = F^,{w**) = -v - F{w*). Since F"{w)w = w* by 
Proposition 1.4 i), we have 

-ly - F^w) ^ F{w*) ^ F{F"{w)w) = F{w) - 2F{w) + 2F{e), 

which is to say that F^,{w) — F{w) — 2F{e) — i^. This implies 

F^{F"{w)x) = F{F"{w)x) - 2F{e) ~v = F{x) - 2F{w) - v. 



where the last equality follows from Equations (|2.6D and (2.7). This concludes the proof 



□ 



Note that together with Equation (2.2), Lemma 2.5 implies that we can replace Axiom (1.5) of the 
original definition of a self-scaled barrier function by the requirement F^,{F" {w)x) = F(x) + C(w) 
for some constant C{w) which depends on w. This fact is already known, see ||l7| . 



3 Group— Theoretic Aspects of Self— Scaledness 

In this section we explore the relationship between the Hessians of self-scaled barrier functions 
and the symmetry group of their domain of definition. Though we present these results primarily 
for the purposes of later sections they are also of independent interest. 



The universal barrier function V(x) defined in ( |l.8| ) plays an important role in the context of this 
section. The choice of the inner product (•, •) used in the definition of the characteristic function 



fxix) via (1.7) is irrelevant, since ipK changes only by an additive constant under a change of 
(•, •). Giiler |^ shows that the universal barrier function U{x) is self-scaled, see Equation (13) and 
Theorem 4.4 in For aU x e int{K) let 

Q{x) ■.^U"{x)-\ 

and let / e int(-fi') be the point characterised by the equation 

Q(/) = i. 

Remark 3.1. It follows from Theorem |l.2| that / is unique. The existence of such a point is also 
well known, see for example page 17 of 

The point / is the "unit" associated with the self-scaled barrier U{x), see Q Proposition 1.3.5, 
p. 14, and it is also the unit of the Jordan algebra associated with U. 

The following lemma is Theorem 3.17, pp. 205-206 in We include a short proof of this 

result because these ideas play an important role in later sections. See also 0, Proposition 1.4.3, 
p. 18 for a different approach to proving this result. 
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Lemma 3.2. The orthogonal subgroup 0(Aut(_ft')) C Aut(_ft') coincides with the stabiliser group 
at f , that is, 

0(Aut(X)) ^{H e AutiK) : Hf = /}. 

Proof. If Ae Aut{K), then 

D^U{Af)[Ah,Ah] ^ D'^Uif)[h,h] 

for every vector he E. That is, A*Q{Afy^A = I, or Q{Af) = AA* , see for example @|, Equation 
(11). This shows that A is orthogonal if and only if I = Q{Af). The uniqueness of / implies that 
this condition is equivalent to Af = f. □ 

Next we note that the elements of Aut{K) have a unique polar decomposition, see Theorem 
3.18, p. 206. For the sake of completeness we give a short proof. 

Lemma 3.3. Let A G Aut(if). There exists a unique vector u G mt{K) and a unique orthogonal 
cone automorphism H G 0(Aut(_ftr)) such that 

A = Q{u)H. 



Proof. By virtue of Theorem 1.2, there exists a unique point u G int(_ftr) such that Q{u)f = 



Af. Then H := Q{u) ^A satisfies Hf = Q{u) ^Af = /, which implies that H is orthogonal 



by Lemma 3.2. Since H is orthogonal and Q{u) is symmetric, A — Q(u)H is indeed a polar 



decomposition of A. 

Suppose now that A = QiHi = Q2H2 where Qi is symmetric and Hi is orthogonal, i = 1,2. 
Then, H :— Q^^Qi — H2H^^ is orthogonal, and we have I = HH* = Q2^Q\Q2^, or = Q\. 
Since Qi and Q2 are symmetric, we have Qi = Q2 and Hi ^ H2. □ 

The following result will play a key role in Section ^ where we classify self-scaled barriers. 
Lemma 3.4. The sets of inverse Hessians of F and U coincide, that is, 

{P{v) : V G int{K)} = {Q{u) : u G int(i^)}, 
and for all x G iiit{K) it is true that 

Fix) = QiQixy/^e~') = Q{xy/^Q{e)-'Q{xy/^ (3.1) 
where G int(i4') is characterised by the equation Q{e^^) = Q{e)~'^. 

The point e~^ is the inverse of e in the Jordan algebra associated with U{x). Note that Propo- 



sition 1.4 Hi) shows that e ^ = —F'{e). 



Proof. If w G mt{K), then Lemma 3.3 implies that we can write P{v) = Q{u)H for some u G 
int(ii:) and H G 0(Aut(ii:)). By the uniqueness of the polar decompositions P(v) — F{v) I and 
P{v) = Q{u)H it must be true that F{v) = Q{u). Thus, 

{P{v) : w G mt{K)} C {Q{u) : u G uit{K)}. 



Conversely, let u G int{K). By Theorem 1.2, there exists a point v G 'mtj K) such that F{v)f 



Q{u)f. But this implies that Q{u)~^P{v)f = f. By virtue of Lemma |3J := Q{u^^)P{v 



is therefore orthogonal. This means that P{v) has the polar decompositions P(v) — F{v) I and 
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P{v) — Q{u)H . The uniqueness part of Lemma 3.2 then imphes that Q{u) = P{v) and H = I. 
This proves the first statement of the lemma. 



Now let X e int(iir) and define u by a; = Q{u)f, see Theorem 1.2. We have 
Qix) = Q{Q{u)f) = Q{u)Q{f)Q{u) = Q{u)\ 



where the second equality follows from the fundamental Formula ( |2.4| ). In a similar vein, taking 
the first part of this lemma into account we obtain 

P{x)^P{Q{u)f)^Q{u)P{f)Q{u). 

These two equations imply that Q{u) — and 

P{x) = Q{xf/^P{f)Q{xf'\ 

In particular, setting x = e yields I = Q{eY''^P{f)Q{eY''^ , that is, P{f)Q{e) = I, and P{f) = 
Q{e)~^ = Q{e^^). The lemma follows, since this implies that 

P{x) = g(.T)i/2g(e-i)Q(x)i/2 = Q(g(x)i/2e-i). 

□ 



Although it does not have a direct bearing on later results, the following proposition already 
shows that the self-scaled barrier function F is intimately connected to the universal barrier 
function. 

Proposition 3.5. There exist constants ai > and ai such that 

U{x) = ai lndetF"(x) + a2. 



Proof. ^From Equation (3.1) we see that det P{x) = det Q(e)~^ det Q{x), implying that det F"{x) = 
ci det U"{x) for some constant ci > 0. Theorem 4.4 in y shows that the function u{x) = hnpK{x) 
satis fies the equation u{x) — C2 + ^ In det u" (a;) for some constant C2. These facts combined with 
(1.8) imply the proposition. □ 



4 Decomposition of Cones and Barrier Functions 



In this section, we prove two related results. A cone is called pointed if it does not contain 
any whole lines. First, we show that any pointed, convex cone decomposes into a direct sum of 
indecomposable or irreducible components in a unique fashion. This result, which is of independent 
interest, is essentially a special case of Corollary 1 in Q . Gruber's paper is the earliest mentioning of 
this result we could find, though it may have been derived several times independently. Gruber's 
original result addresses a more general affine setting which renders his proof more technically 
involved. Therefore, we include a simple and accessible proof. Second, we use this decomposition 
to write any self-scaled barrier function defined on the interior of a symmetric cone as a direct 
sum of self-scaled barriers defined on the irreducible components of K . 

Recall that the Minkowski sum of a set of subsets of E is defined as 

k 

^1 + •• • + \^Xi : Xi G 

i=l 
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If all of the Ai are linear subspaces {0} ^ Ei <Z E which satisfy E — Ei + ■ ■ ■ + Ek and Ei n 
Ej) = {0}, then we say that the sum E = Ei + ■ ■ ■ + Em is direct and write 

E^Ei®E2®---®Era. 



Definition 4.1. Let K C E he a, pointed, convex cone. K is called decomposable if there exist 
cones TO > 2, such that K = Ki + ■ ■ ■ + Km, where each Ki (z = 1, . . . ,to) lies in 

a linear subspace Ei C E, and where the spaces decompose E into a direct sum E — 

-El © -E2 © • ■ ■ © Em- Each iiTi is called a direct summand of X, and X is called the direct sum of 
the {i^i}. We write 

K = Ki®K2®---®Km (4.1) 

to denote this relationship between K and {Ki]'^i. K is called indecomposable or irreducible if it 
cannot be decomposed into a non-trivial direct sum. 



Let us define Ei := (Bj^iEj and Ki (Bj^iKj. If X is the direct sum (O) then every x ^ K 
has a unique representation x — xi + - ■ ■+Xm with Xi €z Ki <Z Ei. Thus, Xi = iTEiX, where tte^ is the 
projection of E onto Ei along Ei. Also, since G Ki,we have = ^i + Sj^^ilO} ^ X^Jli = 
Therefore, 

TTE,K = X, C A'. 
This implies that Ki = tt^.X is a convex cone. Similarly, we have 

(l~nE,)K ^TT^K ^ K,CK. 
We first prove a useful technical result: 

Lemma 4.2. Let K be a pointed, convex cone which decomposes into the direct sum If 
X G Ki is a sum x = Xi + ■ ■ ■ + x^ of elements Xj e K , then each Xj G Ki. 



Proof. We have = ttj, x = Tr^.a^i + . . . + Xk. Each term € Ki C K , therefore we 

have Xj e K and —Xj = X^i^^j ^ Since iiT contains no lines it must be true that Xj = 0, that 
is, Xj = TTEiXj e Ki, (j — 1, . . . , k). □ 



Theorem 4.3. Let K <Z E be a decomposable, pointed, convex cone. The irreducible decomposi- 
tions of K are identical modulo indexing, that is, the set of cones {KiY^Li is unique. Moreover, 
the .subspaces Ei corresponding to the non-zero cones Ki are also unique. In particular, if K is 
solid, then all the cones Ki are non-zero and the .subspaces Ei are unique. 

Proof. Suppose that K admits two irreducible decompositions 

m m 

K = ^K,^^E, and 

i=l i=l 

Note that each non-zero summand in either decomposition of K must lie in span(iir) and that 
the subspace corresponding to each zero summand must be one-dimensional, for otherwise the 
summand would be decomposable. This implies that the number of zero summands in both 
decompositions is codim(span(iir)) . We may thus concentrate our efforts on span(i4'), that is, we 
can assume that K is solid and all the summands of both decompositions of K are non-zero. By 
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Jj), each X £ Cj C K has a unique representation x = xi + ■ ■ ■ + Xm where Xi = iTEiX € Ki C K. 
Also, Lemma 4.2 imphes that Xi G Cj, and hence Xi £ Ki O Cj. Consequently, every x G Cj 



lies in the set (Ki H Cj) + ■ ■ ■ + {Km H Cj). Conversely, we have Ki n Cj C Cj, implying that 
{Kir]Cj)-\ h {Km n Cj) C Cj. Therefore, we have 

Cj = {KinCj) + --- + {KmnC,). 

We have Ki fl Cj C n Fj, Fj = {Ei H Fj) + ■ ■ ■ + {Em n Fj), and the intersection of any two 
distinct summands in the last sum is the trivial subspace {0}. The above decompositions of Fj 
and Cj are thus direct sums. Since Cj is indecomposable, exactly one of the summands in the 
decomposition of Cj is non-trivial. Thus, Cj — KiD Cj, and hence Cj C Ki for some i. Arguing 
symmetrically, we also have Ki C C; for some I, implying that Cj C C;. Therefore, j = I for else 
Cj C Fj O Fl = {0}, contradicting our assumption above. This shows that Cj — Ki. The theorem 
is proved by repeating the above arguments for the cone Ki — (Sk=iiKk = (Bi^jCi. □ 



Next, we show that self-scaled barrier functions have irreducible decompositions as wel l. Let 
F be defined on mt{K) where K is a. symmetric cone with irreducible decomposition ([4.l| ). For 
i — 1, . . . , m, let be a function defined on Ti{Ki), the relative interior of Ki in E. If F{x) = 
S"=i Fi{xi) for every x = X^i^i S ©j^i ri(^i) — int(if), then we say that F is the direct sum 
of the Fi and write 

m 

F = ®F.^- (4.2) 

i=l 



Theorem 4.4. Let K he a symmetric cone with irreducible decomposition (LI). The irreducible 
components Ki are symmetric cones. Let F{x) be a self-scaled harrier for K . There exist self- 
scaled harrier functions Fi for the cones Ki such that 



F ^ Fl' 



>F„ 



The functions Fi are unique up to additive constants. 



Proof. Recall that the universal barrier function U{x) is given in (l.S). Changing the inner prod- 
uct used in the definition of the characteristic function lpk{x) changes U{x) only by an additive 
constant, hence we may assume that {x,y) = '^^i{xi,yi) Ei for the purposes of this definition. 
Here, {■,-)Ei is an inner product defined on Ei chosen so that Ui{fi) = idE^ for some elements 
fi G T:i{Ki) where Ui denotes the universal barrier function defined on vi{Ki). Then we have 
Q{f) — I for / = ©™i/i € i^i{K), in full consistency with our previous notation. Moreover, K 
is self-dual under (■, ■), since K* = Q{f)K = K. Hence, we may choose the vector e S int{K) 
specified in Remark 2.2 as the unique element in iBt{K) such that F"{e) = I under (•,•), see 



Remark B.l 



Thus, U can be written as the direct sum U{x) — (B'^iUi{xi) and Q{x) has block structure 
corresponding to the subspaces Ei, Q{x) = (Bi^iQi{xi) where Qi{xi) = U^^Xi)"^. Consequently, 



Equation (3.1) implies that P{x) also has the same block structure, 

P{x)^Pl{xi)®---®Pm{Xm). (4.3) 

where Pi{x) = Qi{Qi{xi)^/^e~^) £ Aut{Ki) with e~^ = vr^;. (— i^'(e)) , Tr^. being the projection at 
the beginning of this section. 



So far we know that P{x) has block structure corresponding to the direct sum E — ffi™ i-E;, but 
it is not a priori clear that Pi{xi)~^ is the Hessian of a function defined on i-i{Ki). Let the spaces 
Ei be defined as earlier in this section, and let us consider the vector fields Vi : x t—^ ■nEiF'{x), 
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defined on mt{K) and taking values in Ei for all z, i = 1, . . . , m. We claim that Vi depends only 
on Xi = TT^.x. In fact, for any two vectors x,y int(iir) such that Xi — yi we have 



Vt{y) = TTE,F'{y) 



F'{x)+ / F" (ty + {1 ~ t)x)[y - x]dt 



^Vi{x)+ / Pi{TrE,[ty + {I ~ t)x]) TTEily - x]dt ^ Vi{x) + / Orfi, 
Jo Jo 

which shows our claim. Hence, the quotient vector fields 

V, : mt{K)/E, ^ E„ 

X mod Ei i-^ TrEiF'{x) 

are well defined and can be identified with vector fields Vi defined on the cones n{Ki). The direct 
sum of these vector fields amounts to the gradient field 



F' = V, 



int(A') ~> E. 



(4.4) 



F' being conservative, the Vi must be conservative too, implying that these are the gradient vector 
fields of some functions Fi defined on vi{Ki) which are uniquely determined up to additive constants. 
We may choose these constants so that F = (B'^iFi. Clearly, we have Fl'{xi) = Pi{x)~^ for any 
X e mi{K). 



Using Equation (4.3), it is straightforward to check that the Fi are self-concordant, see [|13|. 
Applying Proposition 1.4 i) and iv) to F, using Equation (4.4) and considering variations of 
X S int(ii') only in the part Xi — tt^.x, we obtain that (a;^, —F.[{xi)) = Vi for some number Ui > 0. 
Moreover, applying Proposition 1.4 v) to F and using Equation (4.4) we get Fl{TXi) = T~^Fl{xi) 
for all T > 0. Hence, 



FiiTXi) = 



—F,{^Xi)d^ = F,{x,) + {x^,Fl{^x^))d^ 
F,{xi)- C\x,,-Fl{x,))d(,^F,{x,)-y^ C'd^ 



= Fi{xi) - Vi Inr. 

This shows that the functions Fi are i/i-logarithmically homogeneous. It is a well-known fact 
that any logarithmically homogeneous self-concordant function is also a barrier function, see for 
example [ p^ or |p7| . It remains to show that the functions Fi are self-scaled. As prev iously noted, 
Condition (2^) is satisfied, since Piixi) E Aut(i^i) for all i. Finally, Condition ( |2.6| ) holds for Fi 
because we can apply this condition to F, choosing w — Wi® {(Dj^iei) and x — Xi® {®j^iei) and 
using the block structures of F and F" . 

Note that the irreducible components Ki of K must be symmetric cones, since the Fi are self- 
scaled barriers defined on Ti{Ki). The symmetry of the Ki can also be directly derived from the 
block structure of Q{x) — ®^iQi{xi) and the fact that the set of cone automorphisms {Q{x) : 
X G \nt{K)^ acts transitively on K. □ 



The decomposition Theorem 4.4 shows that for the purposes of classifying self-scaled barriers 
we may concentrate our efforts on irreducible cones. 



5 Classification of Self— Scaled Barriers 



In this section, we give a complete classification of self-scaled barrier functions on the symmetric 
cone K. 
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The definition of a self-scaled barrier function F requires that F changes only by an additive 
cons tant under the action of symmetric cone automorphisms {P{u) : u € 'mt{Ky\, see Equation 
(2.3). However, it is not a priori known how F behaves under the action of an arbitrary element of 
Aut(i^). Note that this is in marked contrast to the case where F is the universal barrier function 
U , which is known to change only by an additive constant under the action of any element of the 
symmetry group of K . This explains in a sense the main difhculty one faces when proving the 
results below. 



The next result is key in resolving this difficulty and is just a slight reformulation of the con- 
jecture raised in [||, according to which self-scaled barriers on irreducible symmetric cones are 
isotropic. Let us denote by Aut(iir)o the connected component of the identity in A\\t{K). 

Theorem 5.1. Let K he a symmetric cone. If H ^ Aut(_ft')o is orthogonal, then F{Hx) = F{x) 
for all X £ int(_fir). 



Proof. Koecher |l^ proves that if if is a symmetric cone, then Aut(if)o is generated by {Q{u) : 
u € V} where V is a neighbourhood of the identity /, see |T^, Theorem 4.9 (b), pp. 88-89. 
Koecher's proof exploits the fact that all derivations of the Jordan algebra associated with U{x) 
are inner. An accessible proof for the case where K is irreducible is given in Lemma VI.1.2, pp. 
101-102, based on certain non-trivial results from the theory of Jordan algebras, li H E Aut(iir)o 
is orthogonal, it follows from Koecher's result that 



H = l[Qiu,) = l[Piv,), 



for some Ui,Vi S mt{K), {i = 1, . . . , Z). Here the second equality follows from Lemma 3.4. There- 
fore, it follows from Equation (2.3) that 



F{Hx) - F{\[ P{v,)x) = F{x) - 2 ^ F{v,) 



Since H f ~ /, setting x — f above yields F{vi) — and settles the claim of the theorem. □ 



The group Aut(i4r)o already acts transitively on K, see 0], page 5. Thus, the above result is 
significant. 

Remark 5.2. Hauser's approach |^ to solving the isotropy conjecture for the cone of positive 
semidefinite symmetric matrices is based on similar ideas. Hauser essentially rediscovers Koecher's 
Theorem 4.9 (b) in this special case, see Proposition 3.3 and Corollary 4.3. He then uses 
Proposition 3.3 in conjunction with the fundamental formula as the key mechanism in the proof 
of the conjecture. Y. Lim |ll[ generalised this approach to arbitrary irreducible symmetric cones 
while refereeing the paper |9| , thus completing the classification of self-scaled barriers. However, 
since Lim was an anonymous referee, his result was not publicly announced. 



Now, assume that K is irreducible and let U be the universal barrier function for K. Let k be 
the rank of the Jordan algebra associated with U{x). Let x be an arbitrary point in 'm.t{K). Then 
there exists an orthogonal frame {ei, . . . , e^} such that x — X]i=i -^jCi, Xi > 0, i — 1, . . . ,k, see 
||l|, Theorem HI. 1.2, pp. 44-45. By C let us denote the cone generated by this frame, that is, 

k 

C={Y^Ke, : A, >0}. 
Note that C is a direct sum of the half-lines {XiCi : A; > 0} and thus a symmetric cone. 
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Lemma 5.3. If F is a self-scaled barrier function defined on the interior of the irreducible sym- 
metric cone K , and if C is as defined above, then 

k k 

F(^aiei) = --'^loga^, ai > (i = 1, . . . , A;). 

i=l i=l 

Proof. Let a be any permutation of {1, . . . , k\. Theorem IV. 2. 5 in implies that there exists an 
orthogonal automorphism H S Aut{K)i^ such that Hci = e„(^i) [i — \, . . . ,k). Using Theorem 5.1, 
we then obtain 

k k 

FC^aiec(i)) ^ FC^aiCi), Va^ > {i = l,...,k). 

i=l i=l 

Define g{ai,... ,ak) '■= F (^'i^i c^iSi) ■ Note that g is a symmetric function. Consider a point 

e + y]'i=i Pi^i ~ T^\=i '^i^i ^ i'i(C'), with arbitrary Pi > and ai :~ 1 + /3i. Applying Theorem 5.1 
in repeatedly, we obtain 



F(e + ^ Ae,) - F{e) = ^(F(e + fte,) - F(e)). 

1=1 i=l 

Using the symmetry of g and F{e) = 0, the above equation translates into 

k 

g(ai,... =^g(aj,l,... ,1), Va, > 1 (i = 1, . . . , A:). 

1=1 

li ai = . . . — ak ^ a above, we have g{a, ... , a) = F{ae) — F{e) — v log a — —v log a. This gives 
.g(a, 1, . . . , 1) = — ^ log a for all a > 1. Consequently, we have 

k 

5(ai, . . . ,Q;fe) = --^logaj, Va^ > 1 (i = 1, . . . , A;). (5.1) 



i=l 



Now, if cki > (z = 1, . . . , fc) are arbitrary, choose t > such that OLi = Oiijt > 1. Since F is 
logarithmically homogeneous, we have g{ai, . . . , ak) = .<?(<ii, . . . , (3;^) — I'logi by the logarithmic 



homogeneity of F. A simple calculation shows that (5.1) holds true for all ai > 0. The lemma is 



proved. □ 

The following theorem classifies self-scaled barrier functions for irreducible symmetric cones. 

Theorem 5.4. Let K be an irreducible symmetric cone, and let F be a self-scaled barrier function 
defined on int(if). Then there exist constants a > and P such that 

F{x) = aU{x) + P, 

where U{x) is the universal barrier function on inti^T. 



Proof. Lemma 5.3 describes the restriction of F on ri(C). Since the universal barrier function is 
also self-scaled, the same considerations apply to U{x) — U{e). Thus, the functions F and U are 
homothctic transformations of each other on ri(C), that is, there exist a > 0, P such that 

F{x) = aU{x) + p. (5.2) 

Let y £ mt{K) be an arbitrary point with the spectral decomposition y = X]f=i '^idi- Corollary 
IV. 2. 7 in ||l[ implies that there exists A e 0(Aut(if )o) such that Aei = di (i = 1, . . . , fc). We have 
y ~ Ax where x — y^,_i VjCi £ ri(C). Theorem 5T gives F{y) — F{x) and U{y) = U{x), and 



hence the Identity (x2) extends to all of \ni{K). □ 
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We are now ready to give the final classification theorem for self-scaled barrier functions on 
arbitrary symmetric cones. This theorem shows that all self-scaled barrier functions are related 
to the standard logarithmic or the universal barrier via homothetic transformations. 

Theorem 5.5. Let F be a self-scaled barrier function for a symmetric cone K with irreducible 



decomposition (4.1). Then there exist constants cq and ci > 1, . . . , Cm > 1 such that 

m 

F ^ cq ~ ^Ciln detKi , 



see Equation ( [4.2[ ). Here det^i Xi denotes the determinant of Xi £ Ti{Ki) in the Jordan algebraic 
sense, see Chapter 2. Conversely, any function of this form is a self-scaled barrier for K. 



Proof. Theorems 4.4 and 5.4 imply that there exist constants do and di > 0, . . . , dm > such that 

F{x) = da + diUi{xi) + . . . + dmUm{Xm), 

where Ui{xi) ~ InipKiixi). It is known that Ui{xi) = const — ni/r.; In det , where is the rank 
of the Jordan algebra associated with Ui{x), and Ui is the dimension of the cone Ki, see 0, 
Proposition III. 4. 3, p. 53. Finally, Theorem 4.1 in implies that the function — alndeta;^ is 
self-concordant if and only if a > 1 . □ 
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